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In brane world scenarios inflationary vacuum energy may escape into the higher dimensional bulk 
leaving behind a dark radiation effect on the brane. The paper analyses the damping of an inflaton 
by massless bulk scalar radiation and the production of dark radiation in a Randall-Sundrum type 
of model. 
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I. INTRODUCTION 

^ : 

In brane world scenarios, matter fields are confined to a four dimensional dimensional hypersurface, or brane, 
| embedded in a higher dimensional bulk. Brane cosmology can reproduce the most important features of the the more 
usual four dimensional cosmology, but sometimes the extra dimensions make themselves known and one such effect 
is considered here. 

Many brane world scenarios are built around the Randall-Sundrum models P, Q, where there is a five dimensional 
anti-de Sitter bulk spacetime and either a double or a single brane. Inflationary solutions to the single brane Randall- 
Sundrum model can be obtained by simply 'detuning' the vacuum energy of the brane 0, • This can be done with 
the extra vacuum energy of a brane or a bulk inflaton field. The former case is called matter field inflation p| . 

In 4-dimensional models, hot big bang radiation is produced by the decay of the inflaton field. In 5 dimensions, the 
inflaton can decay into ordinary matter fields or into bulk fields. In most models of inflation, the coupling between 
the inflaton and other matter fields is very small, and may be smaller even than the coupling to gravity or other bulk 
, fields. Some of the inflaton's vacuum energy can easily end up in the bulk. 
y—{ ' In the most extreme case, nearly all of the vacuum energy of the inflaton field may decay into bulk radiation. This 
7-H , happens according to one scenario, put forward by Enqvist et al. |g,|2|. I n their scenario, the entropy of the universe 
arises from fields associated with the flat directions of the Minimally Supersymmetric Standard Model (MSSM). For 
this scenario to work, the energy released onto the brane by these fields has to dominate over the energy produced by 
■ the decay of the inflaton, which can happen if the inflaton decays into bulk degrees of freedom. 
Qh| The transfer of energy from a brane to the bulk is always associated with a geom etrical effect on the brane which 
JL . mimics a radiation source term in the cosmological evolution equations 0, Ifll iTcl ITU Ilii. IT3. ITEj . The size of 
this dark radiation term in relation to the real radiation is an important feature, because we have strong evidence 
that the real radiation provides the dominant proportion of the radiation density driving the expansion rate at the 
L* ■ nucleosynthesis era. 

The purpose of this paper is to construct a simple model of the inflaton decay process with a single brane and a 
massless scalar bulk field. The bulk field might occur as a hypothetical superpartner of the graviton or a bulk dilaton 
field, for example. We shall construct the reheating terms in the inflaton's equation of motion and solve for the case 
where bulk radiation is the predominant form of reheating. We shall then examine the cosmological evolution to 
determine the residual amount of dark radiation. 
The action for the model takes the form 

S g = -\f (R-2A)dv + -^[ Kdv+ [ £ b dv+ f (C m -a)dv (1) 
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where we have a curved space M. with scalar curvature R and a boundary dM. with extrinsic curvature K. The 
covariant volume integral is denoted by dv. We will take coordinates x M , fi = 0, 1,2,3, along the brane and x 5 = z 
perpendicular to the brane. In the final state after inflation, the vacuum energy a and the 5-dimensional cosmological 
constant A satisfy the fine tuning relationship 6 A = n\u 2 1, 2J. The effective four dimensional Newton's constant G 
is given by 48irG = n\o. 

The inflaton field cj>(x) is restricted to the brane, with Lagrangian density 

C m = ~(d^) 2 -V(<P). (2) 
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In the bulk, there is a massless conformally coupled scalar field with Lagrangian density 

A = -\{da3>) 2 ^m 2 . (3) 

Conformally invariant boundary conditions fix the normal derivative of the Bulk scalar. We consider an interaction 
which allows <fi to spontaneously decay into $. Since the <f> field is confined on the brane, the two fields can only 
interact at z — 0. A simple model for the Lagrangian density of the interaction is 

Cj = -^(xMxf5(z), (4) 

where A is the coupling constant 0. 

The effect of the radiation on the inflaton field equation is a damping term. The general form of this term is 
identical to the four dimenisonal case which can be found in ref. For a flat brane, 

d 2 4> dv 

^+^r J (5) 

where the leading order contribution to the dissipation is 

J(t) = A 2 Im [ d i x'G F (x,x') 2 (j)(t')9(t~t') , (6) 



and Gf{x,x') is the Feynman propagator of the bulk field $ with Robin boundary conditions on the brane. 
For flat extra dimensions, 

Gf(x, x) = /-^-^4cos qz cos qz'^-^^e-^-^, (7) 

where u> = (k 2 +q 2 ) 1 ^ 2 . The intergral © is divergent, but we can apply dimensional regularisation. In five dimensions 
we find 

1 d^ ' - 2 



J V = -^ l »{-,Jt)[jt) +®> (8) 

where fi is the one-loop renormalisation scale. The same renormalisation scale appears in the one-loop corrections 
to the potential. In order for these corrections to remain relatively small, [L has to be large compared to other mass 
scales in potential. 

The logarithmic time derivative appearing in eq. © is a feature of radiation back reaction effects in odd spacetimc 
dimensions [17|. For practical applications, we can use an integral kernel to define the logarithmic derivative (see 

0). 

where 7 is Euler's constant. 

This result which we have found applies to flat branes and flat extra dimensions. However, we can make some 
simple generalisations. In the first place, for the special case of the conformal scalar field, we can use the same result 
for any bulk spacetime which is conformal to flat space. In particular, the result also applies to anti-de Sitter extra 
dimensions. 

Another generalisation can be made if the brane is spatially homogeneous with expansion rate H in the range 
H 2 <C 4>. In four dimensions, the amplitude of density perturbations is of order H 2 /cb, and H 2 <C (f> is consistent with 
having small amplitude density fluctuations. In this case the time evolution of the inflaton is the dominant effect 
producing the dissipative term in the inflaton equation and we can use the flat space result |JSJ. 

The complete form for inflaton equation of motion on an expanding brane with H 2 <C ^ is 

^ +3H ^ + ^ = - cln {^di)^ (10) 

where c = A 2 /(5127r 2 ). The form of the damping term is strongly dependent on the fact that we have not introduced 
any explicit mass terms into the bulk fields. When the same calculation is repeated for a bulk field of mass M, 
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there is no damping of the inflaton when <f> < M 2 . This does not, however, rule out the possibility of higher loop or 
non-perturbative damping effects. 

Inflation is associated with a slow roll approximation in which the leading order terms in the inflaton equation have 
the smallest number of time derivatives. The damping term has too many time derivatives to appear at leading order. 
The effects of damping will start to appear after inflation, in the period which is usually associated with reheating. 
The situation is similar to reheating, appart from the fact that, according to our assumptions, the radiation is radiated 
into the bulk. 

To investigate the period following inflation, consider a quadratic form of potential with minimum at <f> = 0, 

V(4>) = l -m 2 <? (11) 

An approximate solution can be constructed by taking 

(f>(t) =A(t)sin6(t) (12) 
where A(t) is a slowly varying function of t and 9 w rat. From the definition (|9"|>. 

\p at ) at z \ P J 

where Si and Ci are sine and cosine integrals. After solving (|10|) for an initial value (f> = fa at t = t^, we obtain 

A(t)«^(^) 3/2 e — rn(t-u)/i (14) 

If we also set the initial condition to coincide with the end of inflation, then fa ~ M p , where M 2 = (87rG) _1 . 

The effect of the inflaton on the expansion of the universe can be obtained from the brane cosmology equations 
0, E| ITfl ITTL Il2l ITsI ] . We shall only consider the case where the intrinsic brane vacuum energy is large compared to 
the density, a <^ p, when 

QH + 12H 2 = 8irG{p - 3p) - 2k 2 5 T£ 5 (15) 

where T^ b is the stress-energy of the bulk radiation. Since we are assuming that the radiation in mostly into the extra 
dimensions, the density p and pressure p are those of the inflaton field, 

P = \<P 2 + \m 2 tf (16) 
p = \fi - \m 2 ^ 2 (17) 

The radiation term can be obtained by similar means to the dissipative term |T^ . The terms up to order A 2 are 

rpr 3A (ld\d^ dcf> (ld\d 2 cj ) 

where c' = 0(X 2 ). The leading term averages out to zero over a period of the oscillating inflaton field. 
The dark radiation Pdark, defined by 

3H 2 = 87rG(p + Pdark ), (19) 

can be obtained by integrating H15|l. We shall express the result in terms of the radiation flux Tq 5 , using the 
conservation of energy, Tq 5 = —\J<\>. If we assume that Pdark — ► as t — > —00, then 

1 r°° k 2 1 r°° 

Pdaxk = -- / a{t'f T£ 5 dt'~-^- H{t')a{t'fn b dt' (20) 
a J -00 8ttG a 4 J_ OQ 

The first term in identical to the amount of radiation which would be produced if the inflaton's energy was decaying 
into brane radiation rather than bulk radiation. The second term turns out to be less important than the first term 
when a < p. 

In the early stages of reheating, the inflaton field energy density dominates the Friedman equation but decays faster 
than the dark radiation density, which will eventually come to dominate. There is some intermediate time t eq at 



4 



which the scalar field energy density is equal to the dark energy density. In the WKB approximation, after averaging 
over the phase, the radiation flux is given by 

To^icAWlnQ (21) 

The inflaton looses energy provided that m < p.. Most of the the dark energy is generated whilst the scalar field 
dominates eq. lfT§|l. In this regime, 

3H 2 « A-KGm 2 A 2 (22) 
We can solve this equation and evaluate the integral to obtain the dark energy at the time t eq , 

^^-^cV^lnQ (23) 

where M 2 = (87rG) _1 . This is also the maximum value reached by the dark radiation density. 

For comparison, the value of the potential at the end of the inflationary era would be of order Vi ~ m 2 M 2 . 
Therefore Pdark/Vi ~ 10 _8 A . A typical inflationary model, which was conventional in every respect appart from the 
reheating, might have Vi ~ 10 16 GeV. There is no particular reason why A should be very small, and a dark radiation 
'temperature' as high as 10 14 GeV is a possibility. 

We have so far said nothing about the generation of real radiation. Let us suppose that this is produced from a 
field which lies along a flat direction in the potential. These fields have a non-vanishing potential Vp as a result of 
supersymmetry breaking and possible non-renormalisable terms If the dark energy density pdark < Vf at t eq , 
then reheating (or more specifically preheating |20|) would give a radiation energy density p r ~ Vp, and the ratio 

E*«± „ 1(T 8 A 4 — (24) 
p r Vf 

If, on the other hand, pdark > Vf and H ^> Vp 1 / 2 , then the field evolves very slowly along the flat direction. One of 
two things can happen. The expansion rate H can fall below Vp 1 ! 2 whilst Pdark > Vf- The field will then oscillate 
about the minimum of the potential and reheat to give a radiation energy density p r < Pdark- The second possibility 
is that pdark falls below Vf whilst H > Vp 11 / 2 . In this case, we can have reheating with p r > pdark and one has to 
analyse the combination of dark reheating with real reheating for the particular model to obtain the ratios of dark 
radiation to real radiation. 
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